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ABS TRACT 


A turbofan engine propulsion system in which large amounts of water 
are injected into the fan discharge duct is investigated with the goal 
of increasing both the thrust and propulsive efficiency while retaining 
the light-weight qualities of a standard turbofan engine. A parametric 
computer analysis is used to examine the effect of several variables, 
including water-to-gas generator air ratio, water injection velocity, 
fan duct pressure loss, and fan duct thermal and dynamic nonequilibriun, 
upon thrust and propulsive efficiency. In addition, the design param- 
eters of fan pressure ratio and fan bypass ratio are examined for their 
optimum values, and optimum operating combinations of water-to-gas ratio 
and water injection velocity are determined. 

A test apparatus is developed for the direct measurement of wall 
friction force in two-phase flows. A computer program is presented to 
reduce experimental data and compare with pressure drop predicted by 


two empirical correlations. 
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I. COMPUTER OPTIMIZATION OF WATER-AUGMENTED TURBOFAN CONCEPT 


A. SOURCE AND BACKGROUND 

Transportation by water-borne craft is one of the oldest and most 
economical means known to man. Although many improvements have been 
made in payload carrying ability, navigational ability and propulsion 
systems over the centuries, the degree of improvement in the speed of 
water craft has been quite small. It has become desirable, from both 
a military and civilian point of view, to find new light-weight pro- 
pulsion systems for water craft which provide significant increases in 
speed capability while retaining the load carrying ability and effi- 
clency of previous sea transportation systems. 

The feasibility of producing a ship or boat capable of sustained 
high speeds with relatively large payloads and reasonable efficiency 
depends largely upon two factors: drag reduction and the availabil- 
ity of a light-weight efficient propulsion system. Drag reduction 
can be achieved by use of hydrofoils, the captured air bubble concept 
or recent boundary layer control schemes including blowing, suction 
and polymer injection. The development of a suitable propulsion sys- 
tem, however, has not be quite so active. 

The turbofan provides relatively high thrust and very low system 
weight compared to more conventional marine propulsion systems. How- 
ever, in the speed range contemplated for marine craft, 100 knots and 
below, the propulsive efficiency of the turbofan is too low to allow 
serious consideration. 

Early methods of increasing the thrust and propulsive efficiency 
in turbojet and turbofan engines centered on injection of water into 
the compressor or exhaust section. This method of water injection 
has been used primarily in aircraft applications for short periods 
only and has several drawbacks for marine applications, among which 
are stiff requirements on water purity due to machinery corrosiveness 
and energy transfer losses caused by water evaporation due to mixing 
with hot gases. A more recent proposal has been to inject large quan- 
tities of water into the fan discharge duct of a turbofan engine, thus 
overcoming the two primary drawbacks mentioned above because the gases 


are relatively cool and the water contacts no moving machinery. Con- 


Wal 


sequently, sea water can be used in large quantities. 

Several investigators have made preliminary studies concerning the 
feasibility of using water-augmented air-jets for water craft propul- 
sion. Muench and Keith [Ree 1] conducted a preliminary parametric 
study of possible propulsive efficiencies and concluded that reason- 
able efficiencies can be realized. Davison and Sadowski Ref. 2] 
applied their analysis to a particular existing turbofan engine design 
with water injected into the fan discharge duct. Their analysis showed 
that although efficiencies were slightly lower than for other more con- 
ventional marine propulsion systems, the extremely low system weight 
more than compensated for the efficiency deficit. Quandt [Ref 3} 
has assembled some recent theoretical and experimental techniques con- 
cerning "the gas-phase-continuous two-phase jet system." | 

Knudson [Ref. 4] followed the work of Davison and Sadowski in 
proposing the injection of large amounts of sea water into the fan dis- 
charge duct of a turbofan engine to achieve increased thrust and pro- 
pulsive efficiency. The general arrangement of the water-augmented 
turbofan engine is presented in Fig. 1. Knudson was the first to con- 
sider the effects of velocity slip ratios and temperature ratios be- 
tween water droplets and air in the fan mixing duct and nozzle. The 
effects of variation of water- to-gas mass flow ratio and water injection 
velocity were also examined. 

The present analysis undertakes a more detailed computer study 
which includes the pump work required to inject the water at velocities 
greater than craft velocity. A wide range of design and operating param- 
eters is examined for the effect on overall thrust and propulsive effi- 
ciency. Some minor programing errors were discovered in Ref. 4 and 


these are corrected in the present analysis. 


B. WATER-AUGMENTED TURBOFAN ANALYSIS 

The analysis which follows is based on that given by Knudson with 
appropriate changes and corrections for the present use. The analysis 
differs from that for a normal "dry" turbofan in the fan mixing duct 
and fan nozzle sections due to the injection of water. The equations 


are developed in a manner similar to their use in the computer program. 


LZ 


i]s Assumptions 


Following general practice in the analysis of turbofan engines 
the air was assumed to behave as a perfect gas with two constant values 


of Specific ‘heat ratio, x 1.4 for the relatively cool sections of the 


) 
system and 1.33 for the high temperature sections of the engine. In the 
burner section of the engine the arithmetical mean of the two values of 
specific heat ratio was used since the air transitioned from the cool to 
the hot state. It was assumed that no heat was transferred except in 
the combustion chamber and that no work was transferred other than in 
the compressor, fan and turbine. The mixing duct area was assumed con- 
stant and it was also assumed that only converging nozzles were used in 
the engine. 

Finally, it was assumed that completely dry air entered the 
system. The external drag of the water inlet scoop was not included 
in the analysis since this factor is dependent on the actual scoop 
design employed. 

2. Gas Generator Diffuser 


The ambient total pressure is determined from 


Xe 
i Veo-i_ py*) %-1 
oo an : ama, (1) 
and the diffuser efficiency is defined by 
Yo = Fra ~ Pao C2) 
Pro — Fro 


Thus 


Pa = fo t Np (Pro - Pro) (3) 
The total temperature remains constant through the diffuser and is 


%-1 2 
Tit = Tro = Taol! +834 Me) (4) 


The total enthalpy is 


Hrs = Cpe (Tra - Tao) + Hrera (5) 


ie 


where a is the cool temperature specific heat and HREPA is the 


enthalpy of the air at the ambient static temperature, tA0° 


3. Compressor 
The ideal (isentropic) temperature at the compressor exit for 
a specified total pressure ratio, P_,/P is 


TZ 
uo - 4 


Treat = Tra( ”/p., ‘ (6) 


If the compressor efficiency, Ve , is defined as the ratio of the 
isentropic compressor work to the actual compressor work, the total 


enthalpy at the compressor exit is 


Coe (Trax ~ Th 
Hep = —ceelTez-Trs) 4. yy, 


i 
We (7) 
and the corresponding total temperature is 

5, - slg Hee oe (8) 

on 

4. Burner 
The specific heat in the burner section is taken as 
Coc t Coh 

Cpg = a (9) 

The burner exit total enthalpy is 
Hrs = Coa (Tr3 - Tr2) +t Hro (10) 


where the total temperature at the burner outlet, is fixed by the 


t 
Toy 
maximum allowable turbine inlet total temperature. The energy equation 


through the burner can be written 
rc Hra + Mey (Heu2t MeHVE) = (Met p,) Hrs qu) 


where m,, is the compressor mass flow rate, me, is the fuel mass flow 


C 
rate, He is the fuel enthalpy prior to entering the burner, 1s is 
the burner efficiency, HVF is the heating value of the fuel and Hing is 


the totaleenthalpy aitehe bumier outlet. Iisit @e assumed that the 


fuel injection temperature is equal to the air temperature at the 


injection point and that the energy needed to vaporize the fuel is 
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included by choosing the lower heating value of the fuel, then the 


fuel-air ratio can be represented by 


FAR = eu — Cen( Tr3- Tre 
ae (12) 
We Yel(HVF)- Con (Tr3- Tre) 


3 Turbine 
The total work output of the turbine must equal the total work 


input to the other components of the engine. 
mz (Hrs = Hr) a 1. (Hr2 - Hry) + Me, (Hyg Hre )4 my, Wp 


where Mp is the turbine mass flow rate, m,, is the fan mass flow rate, 


ma, is the water mass flow rate and Wp is the water injection pump 


work per unit mass given by 


2 2 
Wp = V7 a= Vz ) (13) 


where Vu and Vo are the water injection velocity and the craft 
velocity, respectively, Tr is the water intake system velocity recov- 
ery factor and Vp is the pump efficiency. The total temperature at 


@aew turbine outlet is 


Ts ladda Trt 4) 


If the turbine efficiency is defined as the ratio of actual work output 


to isentropic work output, then 


Hr3 — Hra 


dy | 
T4I T3 1, e, ers) 
and 
seo alin 
¥,-1 
Fra _ Fear = Rs ("Ye ] , (16) 


6. Gas Generator Nozzle 
The nozzle must discharge to atmospheric pressure (unless 


choked) and thus the ideal static temperature at the nozzle exit is 


tof 
Pp 
Tasr 3 lr4 (4/5, : (17) 


be 


Under assumption of constant. specific heats, the nozzle efficiency can 


be represented by 


V) = aS 

N | e_Ly Cis) 
7g AST 

Thus, for a given nozzle efficiency, the static temperature at the 


nNagegle .erpiat , can be found and the nozzle exit velocity can be 


Tas? 
determined from 


Vas = 2 Con(Trs - Tas) (19) 


The exit mach number is given by 


Ne = Vas i Vas (20) 
as Synge 
If the exit mach number exceeds unity, the computations are adjusted 
to produce sonic velocity at the exit with the appropriate exit 
pressure. 
7. Fan Diffuser and Fan 
The fan diffuser analysis is identical to that of the gas gen- 
erator diffuser. Similarly, the fan analysis corresponds to the com- 
pressor analysis. 
8. Fan Mixing Duct Water Injection Analysis 
The analysis of the water injection process in the fan mixing 
duct considers the possibility of variation in the thermal and dynamic 
equilibrium in the fan mixing duct and the fan nozzle, respectively, 
and static pressure changes between the water injection plane (fan 
mixing duct entrance) and the fan mixing duct exit (stations 7 and 8, 
respectively, in Fig. 1). 
The following terms are used throughout the analysis and are 


defimed Here fom convenienge 


Bypass Ratio (BR) = me/rne 
Water to Gas Generator Air Ratio (WGCR) = ™w/me 
Mixture Ratio (MR) = mw/me 
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The mean velocity of the two-phase flow is defined by 
meVY+ my = (me+tm,) V 
and the mean velocity is thus 
/ = V, + MR( Vw) 
1+#MR 
The momentum equation in the mixing duct is 
(me + Mw)dV = -AdP (21) 

If the area occupied by the air is m r/ Pav and the area occupied by 


the water is Paw the total flow area can be written 
: : MR ) 
A FN PaVa Pw Vor /! 
Combining equations (21) and (22) and assuming a constant area mixing 
duct the resulting equation can be integrated yielding 


Vm 1/7, - b 1 MR 
i > Vy + (Mar fas) (3m) (Bota . aan 


In order to conduct an energy analysis through the mixing duct 


it is necessary to consider the effects of vaporization of water. 


p> He 5 and Beg were determined for the computer using 


The values H 
least squares cubic approximations of the data in the steam tables 
Ref. 5 for the temperature range 510-660 degrees Rankine. The 
equations used are shown at the beginning of the computer program in 
Appendix II. The values of -e and - are found from the equations 


=H, + 
H, =H, + He, 


iI 
* 
+ 
op) 


tie OT 


The partial pressure of water vapor, P in atmospheres is 


Pye 
determined from equation (12) in Ref. 5. 


- _ X |A+BY¥ + CX 
logyy (Poy) 2 logyo( Fe ) oO i (24) 
The specific humidity, X, of the mixture is defined as the ratio of the 


mass of the water vapor to the mass of the dry air in the duct. Using 


the perfect gas relations for air and water vapor it can be shown that 


X = _ Feve MWy_\{_Tag (Oss 
8 Faas - Pove MWa Tyg 


ey 


where MW is the molecular weight and subscripts V and A signify water 
vapor and air, respectively. The vapor temperature is assumed equal 
to the air temperature. 


The energy equation through the mixing duct is 


2 
me (Haq? Yea) + Yiwe(Hey # Mga) = ™e (Hag * *) 


2 
(26) 


The analysis thus consisted of the equations (23), (25) and (26) with 


six unknowns: Vig? Vag» Tyg > Tag: Pag and Xe. A solution is found 


by specifying values for three ratios of unknowns: 


TR8WA = Trig/Tag (27) 
VR8WA = Viig/Viag (28) 
PR87. = Pag Pay (29) 


By varying these ratios, the effects of departure from thermal or 

dynamic equilibrium at the mixing duct exit and static pressure change 
through the duct can readily be examined. Although six equations with 
six unknowns exist, an explicit solution is not possible because equa- 


tions (25) and (26) contain involved functions of T through the en- 


thalpy and partial pressure terms. However, it is cml to deter- 
mine Vag and Vig using equations (23), (28) and (29). The remaining 
three equations, (25), (26) and (27), are then manipulated to obtain 

a function of Tag only (called FUN). A Newton-Raphson iterative solu- 
tion is then used in the form 


Taj = Tagj-1 — "UN /DFUN (30) 


ie >. .th , , , a 
where subscript j indicates the j approximation, j-l indicates the 


th 
(j-1) approximation, etc., and DFUN is the derivative of FUN with 


4 


both functions being evaluated at T mee connect 


respect *"to™l A8j-1 


A8? 
Value or Tag is obtained when the absolute value of FUN/DFUN becomes 


zero. It is then a simple matter, knowing Tag» to determine Xo and 


Ti from equations (25) and (27). 
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9. Fan Nozzle Analysis 


Following reasoning similar to that in the preceding section 


the specific humidity at the fan nozzle exit is 


__Pevo v Th 
Z Pv MW A9 
X9 mt) | MW, }( Tyg cx 


The constant area assumption made to allow solution of the 
momentum equation for the fan mixing duct analysis obviously does not 
apply through the nozzle. The approach taken (as in Ref. 2) is that 


the net change in entropy through the nozzle is zero. 


(Sg - Sa) 


Equation (32) when expanded gives 


rg [Coe dm (FE ] _ %)| . (49° Korine) Spot ™eXg Sao] 


Ss.) = (0 (32) 


» (So ~ “8’WATER 


AIR 


~[ (rw Xg™5) Seg + mg Xg S96 va 


The exit pressure, must equal atmospheric pressure, leaving two 


Pag? 


equations, (31) and (33), with unknowns TQ? T49 and Xq- By specifying 


TROWA = Two/r,. e 


and solving equations (31) and (33) for a function of Tag» a Newton- 


Raphson iteration again provides the solution for Tg? Tig and Xo 


The energy equation is 


: . 
mr (Aas z Mas ¢ (rw Xe) (eg t Ys. t+ M, Xp (Hog t 4s 


= rie (Hao ¢ Y42) (tya- Xo7e)(Heg + 8) i Xollgg + 42) (35) 


The fan nozzle exit velocity, Vag: is determined by specifying 


VROWA = Visa! Vag (36) 


and solving equation (35) for Vag as shown on the following page. 


Lg 


l/ - Cec (Tas Tag) if Xa Hege E Xo Heag it MR (Heg- He) 
me [14 Xq + (MR- Xq)( VROWA)*] 






Vis, La (MR- Xo 
[4 + Xg + (Me- Xg 


(vR8wa)* 
\VROWAY 


10. Specific Thrust and Propulsive Efficiiengy 


The thrust of the water-augmented turbofan is determined from 





the momentum flux equation 


TH = m Vas i. Me Vo if (ain ¢ t Xo me) Vag = m, Vo 


te (m1 we -Xq ™me )Vwg- Mw Vwo + A. (Pas - Pro) (38) 


where the water velocity entering the inlet scoop, Vio? is assumed 


equal to the craft velocity, V The specific thrust is the thrust 


9° 
divided by total air mass flow through the engine. 


moet TH - 
ST = mice = walla FAR Vpe- -V, + BR (1+X5 Vig 


14% 
= BRIV.) + (WGR- Xo BR) Vw9 -WGR (Vwo)* Bae Vag As "fp 2) (39) 


The propulsive efficiency is a measure of how well the system converts 
the kinetic energy change of the working medium into thrust power and is 


is represented by the equation 


1 TH (Vo) 
PR TH(Vo) +! + 4m, (Viger V vi + @(m gt Yom. (Vase Vey 4 (wg Xam eX Vg Vwo) eo 


C. OBJECTIVES OP fae PARAMEITRS SLUDY 

The major objective of the present analysis was to refine the com- 
puter procedure initiated by Knudson so that a more comprehensive para- 
metric optimization could be carried out. The primary performance in- 
dicators were taken to be the specific thrust and propulsive efficiency. 
The specific thrust and propulsive efficiency of the water-augmented 
turbofan were compared to their counterparts for a “dry" turbofan which 


provided a thrust ratio and an efficiency ratio for the wet-to-dry case. 


20 


The system parameters described below were varied over ranges of prac- 
tical values to determine effects on thrust and efficiency ratio. 
1. Fan Total Pressure Ratio and Bypass Ratio 
Perhaps the two most important parameters in the design of a 
turbofan engine are the fan total pressure ratio (FPR) and the bypass 
ratio (BR). Consequently it was expected that these parameters would 
also be of importance in the design criteria for a water-augmented 
turbofan engine and that information concerning the effect of bypass 
ratio and fan total pressure ratio related to various water-to-gas 
ratios and various water injection velocities could be very useful in 
optimizing the design of.the fan section. 
2. Water-to-Gas Ratio 
The effect of water-to-gas ratio variation was expected to be 
important in optimizing the thrust and efficiency ratios since the 
basis of the entire concept is water-augmentation. The range of values 
of water volume fraction for mist flow in the fan duct was assumed from 
zero to about 15 percent at standard conditions which corresponds to a 
water-to-gas ratio range of zero to 500. 
3. Water An jeetiongve locity 
For a given water-to-gas ratio a wide range of possible water 


injection velocities (V ) exists. The effect of losses caused by the 


INS 
mixing process in the fan mixing duct is directly related to the veloc- 
ity ratio between the water and air in the injection plane making the 
water injection velocity a significant parameter. The forward velocity 
of a vehicle through the water provides only a certain maximum injection 
velocity without providing additional pump work. If water injection 
velocities are desired which are greater than the maximum ram velocity 
available, allowing for inlet duct losses, pump work must be added 

which necessarily reduces the efficiency and thrust of the gas genera- 
tor. Only if the advantages of increased injection velocity outweigh the 
losses associated with the pump work required will the operation be 
worthwhile. Although the previous work by Knudson allowed variation 

of water injection velocity the pump work required was not considered. 
One of the objectives of the present analysis was to provide the capabil- 
ity for handling injection pump work as required by variations in water 


injection velocity. 
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4. Optimum Combinations of WGR and Ving 


The efficient operation of a water-augmented turbofan propul- 
sion system would require a knowledge of the optumum combinations of 
water-to-gas ratio and water injection velocity to provide maximum 
thrust at a given ship speed. With this in mind it was decided to 
provide information relating the optimum combinations required. 

5. Mixing Duct and Nozzle Thermal and Dynamic Equilibrium 

The water and air enter the mixing duct at different velocities 
and different temperatures, in general. During the mixing process the 
temperatures and velocities, respectively, tend to equalize so that 
different ratios of water temperature to air temperature and water 
velocity to air velocity can occur at the mixing duct exit. Similarly 
in the fan nozzle the accelerations and temperature changes undergone 
by the two phases can result in various temperature ratios and velocity 
ratios at the nozzle exit. It is of interest to determine how these 
ratios affect the thrust and efficiency of the propulsion system. 

6. Fan Duet Pressure Drop 

The change in static pressure of a flow in a duct is directly 
related to the wall shear force through the momentum equation. It is 
possible, therefore, to gain some understanding of the effects of wall 
shear friction losses in two-phase flow by looking at the effects 
caused by variation of the ratio of the static pressure at the fan 
mixing duct exit to the static pressure at the entrance to the fan 


mixing duct. 


D. COMPUTER PROGRAM AND MODIFICATIONS 

The computer program used in the present study is based upon a 
standard turbofan analysis which was modified by Knudson to consider 
the injection of water into the fan discharge duct. An explanation 
of various input and output parameters is given in Appendix I. A 
copy of the program itself and a sample of the computer output are 
presented in Appendix II. 

Originally the computer program consisted of two separate parts: 
a "wet" turbofan program and a "dry'' turbofan. A considerable amount 
of program revision was required to enable the combination of the two 


programs so that the ratios of water-augmented thrust and efficiency 
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to dry thrust and efficiency, respectively, could be calculated direct- 
ly. Of less importance, the output format was modified to present the 
important data in more compact form to allow wider ranges of the var- 
ious parameters to be covered without excessive use of paper. 
ij lotalginmenopy Detenmination at Surnerstumlet 

In the original program the calculation used in determining 
the total entropy at the burner outlet was based on the change in en- 
tropy from the burner entrance to the burner exit added to the refer- 
ence entropy of the atmosphere. However, the use of the reference 
entropy of the atmosphere in Knudson's analysis was incorrect and the 
reference entropy was corrected to be the total entropy at the 
burner entrance. 

2. Temperature [terations 

Several items in the iterations to determine the mixing duct 
air temperature and the fan nozzle air temperature were improved in 
accuracy or were added because of omissions in Knudson's program. 
Of these, the most important was the accuracy indicator used to deter- 
mine the convergence of the Newton-Raphson iteration for air tempera- 
ture at both the duct exit and the nozzle exit. Knudson examined only 
the numerator of the convergence term whereas the proper term to ex- 
amine is the absolute value of the entire convergence term. This was 
corrected for the present analysis and repeatable results were ob- 
tained which in some cases differed significantly from those obtained 


by Knudson. 


E. RESULTS AND DISCUSSION 

The significant results of the computer study are presented in 
graphic form in Figs. 2 through 22. For all results shown the craft 
velocity was held at 50 knots although the computer program can handle 
any speed range. The constant velocity allowed more meaningful com- 
parisons to be made among the results. 

The following parameters were also held constant throughout the 
analysis, the values used being typical of presently available turbo- 
fan engines. The compressor total pressure ratio was maintained at 
13.8 and the turbine inlet total temperature was fixed at 2500 degrees 


Rankine. The dry fan duct total pressure ratio was taken as unity 
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since no separate fan duct exists on a dry turbofan engine. Except 
for the study of duct pressure ratio effects the wet fan duct static 
pressure ratio was arbitrarily fixed at 0.95 since the expected two- 
phase pressure drop was not known. 

The maximum available velocity for water injection without injec- 
tion pump work depends upon duct and piping losses in transferring the 
water from the inlet scoop to the injection nozzles. The velocity 
recovery factor chosen for the present analysis was 0.8. The mach 
number at the fan outlet (water injection plane) was assumed to be 
raz The efficiencies assigned to the various components are listed 


below and are the same as those used in Refs. 2 and 4, 


Diffuser 0795 
Compressor 0790 
Burner 1.00 
Tur bine 0.85 
Nozzle UES 
Fan Diffuser oe 95 
Fan 0.90 
Fan Nozzle O:9e 
Water Injection Pump 0.90 


1. Effects of Fan Pressure Ratio 

It is of interest to determine the design criteria for maximiz- 
ing thrust ratio with respect to fan pressure ratio and bypass ratio. 
The variation of thrust ratio with fan pressure ratio, while holding 
bypass ratio constant at 4.0, is presented for various water injection 
velocities in Fig. 2 and for various water-to-gas ratios in Fig. 3. 
The design envelope shown in both graphs is determined by the amount of 
turbine work required to drive the compressor, fan and water injection 
pump. The turbine work available is dependent on the maximum pressure 
drop allowable across the turbine. As the turbine work required ex- 
ceeds the turbine work available the turbine exit pressure must go be- 
low atmospheric pressure. Since this is impossible, the design limit 
is reached when the turbine exit pressure equals atmospheric pressure. 

It is observed from Fig. 2 that for a given water injection 
velocity or water-to-gas ratio the highest possible fan pressure ratio 


provides the maximum thrust ratio, except for an injection velocity of 
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100 knots, below a fan pressure ratio of 1.1. The increase in thrust 
ratio for low values of fan pressure ratio at high water injection 
velocities and high water-to-gas ratios is considered to be unreal- 
istic because of the magnitudes of certain of the parameters. Specif- 
ically, a fan pressure ratio of 1.0 corresponds to a duct in the free 
stream into which water is injected with no losses occurring in the 
fan. Under this condition higher water injection velocities reduce the 
velocity differential between the water and the faster air thus poss- 
ibly reducing mixing losses sufficiently to allow an increase in 
thrust ratio, 

The optimum fan pressure ratio for duct and nozzle thermal 
and dynamic equilibrium, respectively, appears to be in the range 1.4 
to 1.5 for a design speed of 50 knots. This range gives the maximum 
thrust ratio while still allowing the widest range for both water in- 
jection velocity and water-to-gas ratio. A fan pressure ratio of 1.5 
was used for the remainder of the present analysis. The effect of 
variation of fan pressure ratio on efficiency ratio was also studied 
and exhibited trends generally similar to the corresponding thrust 
ratio curves. | 

? Effects: of Bypass. Ratio 

The selection of the optimum bypass ratio is also an important 
design criteria for the water-augmented turbofan engine. The variation 
of thrust ratio with bypass ratio, while holding fan pressure ratio 
constant at 1.5, is presented in Fig. 4 for various water injection 
velocities and Fig. 5 for various water-to-gas ratios. Again the de- 
sign envelope reflects the maximum available turbine work limited by 
turbine exit pressure and atmospheric pressure. 

Increases in both injection velocity and water-to-gas ratio 
provide increased thrust ratio within the design envelope. Increases 
in bypass ratio generally increase thrust ratio for a given injection 
velocity or water-to-gas ratio except for injection velocities above 
about 70 knots and water-to-gas ratios greater than 400. The decrease 
in thrust ratio with increasing bypass ratio for a given injection 
velocity above about 70 knots (see Fig. 4) may be attributed to the 
fact that the increase in pump work required to maintain the water- 


to gas ratio constant as bypass ratio increases detracts more from 
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the thrust of the gas generator than the amount of thrust increase of 
the fan section. A similar explanation applies to the results of 
Fig. 5 regarding various water-to-gas ratios. 

The optimum design bypass ratio for maximizing thrust ratio 
at a design speed of 50 knots appears to occur in the range of bypass 
ratios from 3.0 to 4.0. This range allows a wide selection of injec- 
tion velocities and water-to-gas ratios while still remaining within 
the design envelope. The bypass ratio chosen for the remainder of the 
present analysis was 4.0. The effect of variation of bypass ratio on 
efficiency ratio followed trends generally similar to the correspond- 
ing’ thrust ratio eurves. 

3. Effects of Water-to-Gas Ratio 

The variation of thrust ratio with water-to-gas ratio is pre- 
sented in Fig. 6 for various values of water injection velocity. 
Figure 7 shows similar information with respect to efficiency ratio. 
It is apparent that both water-to-gas ratio and water injection veloc- 
ity greatly affect the thrust ratio and the efficiency ratio. A maxi- 
mum thrust ratio of slightly less than 4.0 is achieved for a water-to- 
gas ratio in the range 375 to 400 at a water injection velocity slight- 
ly higher than 60 knots and a craft velocity of 50 knots (see Fig. 6). 
However the efficiency ratio shown in Fig. 7 continues to increase 
for values of water-to-gas ratio exceeding 500. At a water-to-gas 
ratio of 500 the maximum efficiency ratio obtainable is about 6.2 for 
a water injection velocity slightly below 60 knots and a craft velocity 
of 50 knots. These results show that. it is indeed advantageous to in- 
crease the injection velocity even at the expense of the pump work 
required, 

The operating envelope shown in Figs. 6 and 7 indicates the 
maximum turbine work available for driving the compressor, fan and 
water injection pump and arises because the turbine exit pressure is 
lower than atmospheric pressure at the higher fan pressure ratios. 

The explanation given earlier in Section D.1 for design envelope ap- 
plies to the operating envelope as well. 

Each water injection velocity has an associated water-to-gas 
ratio at which thrust ratio is maximized and another water-to-gas ratio 


for which efficiency ratio is maximized. This may be explained by the 
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consideration that for a fixed water injection velocity the initial 
velocity ratio between the water droplets and the air is constant. As 
the water-to-gas ratio increases the thrust from the fan nozzle is in- 
creased but the injection pump work increase required causes a decrease 
in the thrust of the gas generator. Initially the thrust increase from 
the fan section predominates but eventually the thrust decrease from 
the gas generator section becomes larger than the thrust gain and the 
net thrust begins to decrease. 

4, Operating Curves 

An important question concerns the optimization of water injec- 
tion velocity and water-to-gas ratio to provide maximum thrust ratio at 
a given craft velocity. A suggested form for the presentation of such 
information is shown in Fig. 8. For a desired craft velocity and water 
injection velocity the water-to-gas ratio required to provide the maxi- 
mum thrust ratio can be determined from the curve. 

The curves do not, however, provide information concerning the 
optimum combination of both injection velocity and water-to-gas ratio 
required to provide the maximum possible thrust ratio at a given craft 
velocity. For example, at a water-to-gas ratio of 200 and craft veloc- 
ity of 50 knots there exist two injection velocities, 41 knots and 79 
knots, respectively, which will provide a maximum thrust ratio, but the 
greater of the two thrust ratios is not specified. However, comparison 
with Fig. 6 indicates thet the higher thrust ratio coincides with the 
higher injection velocity. Figure 7 also shows this case to be the most 
efficient. This suggests that in cases where two possibilities exist 
for injection velocity the higher injection velocity should be used to 
provide the higher thrust and efficiency ratio. 

The importance of the peaks at high water-to-gas ratios on the 
curves for craft velocities of 25 and 50 knots, respectively, lies in 
the fact that, for water-to-gas ratios above about 500, the fan duct 
contains such a high percentage of water that the water-jet principle 
is approached with its increased pumping requirements and consequently 
increased weight. At some point the weight increase will more than 
offset the performance advantage and the water-augmented turbofan loses 
its appeal. For this reason it is suggested that the water-to-gas ratios 


be kept as low as possible. 
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3. Fan Duct and Nozzle Equilibrium Parameters 
One of the more important aspects of the present computer anal- 
ysis has been the determination of the effect of fan duct and nozzle 
equilibrium parameters upon the thrust ratio and the efficiency ratio. 
The important equilibrium parameters are taken to be the fan duct exit 
velocity ratio, V/V the fan duct exit temperature ratio, T/T 
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fan nozzle exit velocity ratio, Vig! Vag» and the fan nozzle exit tem- 


perature ratio, Tg/t these ratios are referred to in the succeeding 


AQ? 
test as, respectively, duct velocity ratio, duct temperature ratio, 
nozzle velocity ratio and nozzle temperature ratio and are water-to- 
air ratios. 

These parameters are important for several reasons. First, 
little is known about the nature of the mixing process; consequently, 
analysis must be made on an empirical basis. Second, Lockhart and 
Martinelli [Ref. 6] and Chenoweth and Martin [Ref. 7] indicate that 
mixing duct pressure losses are greater for the flow of a two-phase 
mixture than for either phase flowing separately. The extent of the 
increase and the factors which govern the increase are not very well 
known. Further discussion of this subject is presented in Section II 
of this paper. 

a. Fan Duct Velocity Ratio 

The variation of thrust ratio and efficiency ratio with 
duct velocity ratio is presented in Figs. 9 and 10, respectively, for 
various values of water injection velocity. Under the restriction that 
all other equilibrium parmeters have the value unity, the duct velocity 
ratio has only a very slight effect on thrust ratio and efficiency ratio. 
The duct velocity ratio, therefore, does not appear to be a significant 
parameter affecting thrust ratio and efficiency ratio. 

b. Fan Duct Temperature Ratio 

The effects of duct temperature ratio on thrust ratio and 
efficiency ratio are shown in Figs. 11 and 12, respectively, for various 
water injection velocities. The initial conclusion from Fig. 1l is that 
it would be highly desirable to maintain the duct temperature ratio as 
low as possible to attain the maximum thrust ratio. However, further 
consideration of the restrictions imposed, namely duct velocity equi- 


librium, nozzle velocity equilibrium and nozzle temperature equilibrium, 
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create unrealistic, if not impossible, conditions for the flow to meet, 

First, it is extremely unlikely that a low value of duct 
temperature ratio could, at the same time, produce equilibrium of the 
other three parameters. Seconds.the experimemial results of Reese and 
Richard [Ref. 8| indicate that thermal equilibrium occurs in two-phase 
flows with contact times as low as the order of one millisecond. A typ- 
ical expected contact time of the order of 20 milliseconds suggests 
thermal equilibrium at all times in the fan duct exit. 

A further consideration is that for low values of duct tem- 
perature ratio the bulk of the mixing process must occur in the nozzle, 
Since the calculations through the nozzle include the assumption of no 
net entropy change, the thrust ratio resulting is larger than that which 
could realistically» occur. These reasons suggest that the results given 
in Figs. 11 and 12 do not give an accurate indication of realistic ef- 
fects. 

c. Fan Nozzle Velocity Ratio 

Figures 13 and 14, respectively, present the effect on 
thrust ratio and efficiency ratio of varying the ratio of water velocity 
to air velocity at the nozzle exit for various water injection veloci- 
ties. For a given injection velocity neither thrust ratio nor efficiency 
ratio varies a great amount, however, the variation is greater than the 
Similarscase for variation of duct velocity ratio. For an injection 
velocity of 100 knots the thrust ratio increases from about 3.13 at a 
velocity ratio of 0.5 to about 3.22 at a velocity ratio of 1.0 and the 
efficiency ratio increases from about 3.77 to about 3.88. Both curves 
become nearly horizontal near a duct velocity ratio of 1.0. 

It should be pointed out that the restriction of equilib- 
rium conditions for the other parameters is not so unrealistic at the 
nozzle exit as it was at the duct exit. This is because the flow is 
accelerated through the nozzle and the air, because of its lower density, 
accelerates more quickly than the water causing a velocity difference at 
the nozzle exit regardless of the equilibrium state at the fan duct exit 
(nozzle entrance). As was the case with duct velocity ratio, the nozzle 
velocity ratio appears to have little significant effect upon thrust 


Babiovand efiiiciency rajgio. 
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d. Fan Nozzle Temperature Ratio 

The variation of thrust ratio versus nozzle temperature 
ratio is presented for various water injection velocities in Fig. 15, 
and for various water-to-gas ratios in Fig. 16. In both cases the 
thrust ratio decreases smoothly with increasing nozzle temperature ratio 
indicating that the most desirable operating range is near thermal 
equilibrium where the curves are nearly flat. Again, based on the re- 
sults of Reese and Richard, it is probable that the nozzle temperature 
ratio will be near unity in spite of the temperature changes of the air 
through the nozzle. 

The variation of efficiency ratio with nozzle temperature 
ratio is shown in Figs. 17 and 18, respectively, for various water in- 
jection velocities and various water-to-gas ratios. The results tend 
to follow trends similar to those for the thrust ratio except for the 
higher values of water injection velocity and water-to-gas ratio. For 
a water injection velocity of about 80 knots in Fig. 1/7 the efficiency 
ratio is relatively unchanged over the range of nozzle temperature 
ratios. An injection velocity of 100 knots shows increasing efficiency 
ratio with increasing nozzle temperature ratio. Similarly a water-to- 
gas ratio of 400 or higher in Fig. 18 exhibits a similar trend of in- 
creasing efficiency ratio. 

The information for temperature ratios. significantly greater 
than unity is considered to be unrealistic in view of the predicted short 
contact times required to achieve thermal equilibrium. In addition the 
large pressure drop through the nozzle, which would be required for such 
a great air temperature change, is highly unlikely to occur. 

6. Effects of Fan Duct Pressure Ratio 

An important parmeter affecting the performance of the water- 
augmented turbofan engine is the pressure loss in the fan mixing duct. 
Since the wall shear force in two-phase flow is an unknown quantity 
which has not been treated analytically in this paper, a range of values 
for the ratio of static pressure at the mixing duct exit to the static 
pressure at the duct entrance has been used to determine the effect on 
thrust ratio and efficiency ratio. Figures 19 and 20 present thrust 
ratio versus duct pressure ratio for various water injection velocities 
and various water-to-gas ratios, respectively. Figures 21 and 22 pre- 


sent parallel results for efficiency ratio, respectively. 
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It is theoretically possible for the static pressure to in- 
crease in the direction of flow if large velocity variations exist over 
the injection plane and the wall shear stress is small. nus, these 
plots have been extended to preSsure ratios greater than one. It is 
easily seen thet the most desirable range of operation to provide max- 
imum thrust ratio is at the maximum duct pressure ratio. This is 
simply the expected result that minimum duct friction provides the 
best operating conditions. Also, as predicted by a momentum analysis, 
the highest water injection velocity provides the least sensitivity to 
duct pressure losses as well as providing the best thrust ratio. 

The effect of duct pressure ratio on efficiency ratio is simi- 
lar to the effect on thrust ratio; however, for a fixed water-to-gas 
ratio the efficiency ratio is virtually unchanged with varying duct 
pressure ratio. Again increasing water injection velocity decreases 
the effect on efficiency ratio of varying duct pressure ratio. As the 
water injection velocity increases to about 80 knots, the injection 
pump work required begins to predominate over the efficiency gain af- 
forded by the high water injection velocity, causing the efficiency 


ratio to decrease with increasing duct pressure ratio (see Fig. 21). 


F, CONCLUSIONS 

The injection of water into the fan discharge duct of a turbofan 
engine provides a feasible propulsion system for high speed sea-borne 
vessels, Among the most important parameters affecting the thrust 
and propulsive efficiency of water-augmented turbofans are the water- 
to-gas ratio and the water injection velocity. Each water injection 
velocity has an associated water-to-gas ratio which maximizes thrust 
and another water-to-gas ratio which maximizes propulsive efficiency. 
thrust ratios and efficiency ratios of the water-augmented turbofan 
to the dry turbofan of three or four are easily attainable. It should 
be noted that these figures do not account for the increased drag that 
would be caused by the water inlet scoop. 

Velocity slip ratio in either the tan mixing duct or the fan nozzle 
has negligible effect upon the thrust ratio and efficiency ratio. The 
effect of temperature difference between water and air is somewhat 
greater than for the velocity ratio. However, the likelihood of non- 


equilibrium for either fan duct or fan nozzle temperature ratio is 
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slight due to the short contact times required for thermal equilibrium. 

Thrust ratio and efficiency ratio both increase with increasing 
fan duct pressure ratio. Fan duct pressure ratio can be related to 
fan duct wall friction forces through the momentum equation, indicating 
that minimum wall friction is desirable. Thrust ratio and efficiency 
ratio become less sensitive to fan duct pressure ratio with increasing 
water injection velocity and increasing water-to-gas ratio. Since 
little is known about two-phase pressure drop, experimental work should 
be done to increase the knowledge in the field. Section II of this 
paper describes the development of a test facility for wall shear force 
and pressure drop measurements in two-phase flow. 

For a design speed of 50 knots, the ideal fan total pressure ratio 
falls within the range 1.4 to 1.5, if use of a wide range of water 
injection velocities and water-to-gas ratios is desired. At the same 
design speed, the optimum fan bypass ratio occurs in the range 3.0 to 
4.0, to allow use of a wide range of water injection velocities and 
water-to-gas ratios. 

For any given water injection velocity a set of curves may be con- 
structed which provide the water-to-gas ratio required to achieve maxi- 
mum thrust for a specified craft velocity. The present analysis has 
not considered the relationship between increasing the water injection 
velocity and the resulting weight increase due to pump requirements. 

As the water injection velocity and water-to-gas ratio increase the 
pure water-jet concept is approached and the performance advantages 
of the water-augmented turbofan concept, primarily the low system 
weight, are lost. A further study is thus required to determine the 
point at which weight and cost increases for water injection pumps 


begin to detract significantly from the overall system effectiveness. 
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II. DEVELOPMENT OF A TEST FACILITY FOR TWO-PHASE FLOW 


A. TWO-PHASE PRESSURE DROP PREDICTION 

thé ability to predict the pr@ssure drop occurring in a two-phase 
duct flow is vital to the accurate analysis of the water-augmented | 
turbofan concept discussed in Section I. To the author's knowledge no 
complete, accurate, theoretical analysis of two-phase pressure drop 
has as yet been proposed. However, Dukler, Wicks and Cleveland [Rez. 9] 
and other available literature indicate that two empirical data correl- 
ations for the prediction of two-phase pressure drop have given con- 
sistently better results than other correlations. These correlations 
were developed, respectively, by Lockhart and Martinelli [Ref. 6 and 
Chenoweth and Martin [Ref. 7| . An explanation of each of the correl- 
ations and some of the special assumptions made in them is presented 
in the following text. 

1. Lockhart-Martinelli Two-Phase Flow Correlation 

The correlation of Lockhart and Martinelli is dependent upon 

which of four types of flow exists during the simultaneous flow of a 
liquid and a gas (or vapor). These flow regimes are: 


1. Flow of both the liquid and the gas may be turbulent 
(turbulent-turbulent flow) 


2. Flow of the liquid may be viscous and flow of the gas 
may be turbulent (viscous-turbulient flow) 


3. Flow of the liquid may be turbulent and flow of the gas 
may be viscous (turbulent-viscous flow) 


4. Flow of both the liquid and the gas may be viscous 
(viscous-viscous flow) 


"Viscous" flow is the term used by Lockhart and Martinelli for what 
is more properly called laminar flow. 

The assumptions made are that the "static pressure drop for 
the liquid phase must equal the static pressure drop for the gaseous 
phase regardless of the flow pattern, as long as an appreciable radial 
static pressure difference does not exist' and ''the volume occupied 
by the liquid plus the volume occupied by the gas at any instant must 
equal the total volume of the pipe." It is further assumed that the 
transition from laminar to turbulent flow occurs in the range of 


Reynolds numbers between 1000 and 2000. The Reynolds number is calcu- 
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lated using the mass flow rate of the phase in question and the inside 
diameter of the pipe. For example, the liquid phase Reynolds number is 


determined from the equation 


B 
Re, = a = (41) 


where ma is the liquid mass flow rate, in pounds-mass per second, 
D is the inside pipe diameter in feet and M, is the absolute viscosity 
of the liquid in pounds-mass per foot-second. 

Having determined which of the flow regimes is appropriate 
from the Reynolds number calculation for each phase, the all-liquid 
and all-gas pressure drops, respectively, are calculated. The all- 
liquid pressure drop in a pipe L feet in length is based upon the 


liquid flowing at a rate ma with a density (, as in the equation 


oy 
AFP, = PL Lo (42) 


where the liquid velocity, V is determined from continuity, assuming 


oe 


that the liquid occupies the entire pipe, and f, is the liquid friction 


L 
factor based on Re, . Similarly, an all-gas pressure drop is found from 
2 
Be Lo VG 
Wee - Be ba (43) 


The ratio of the all-liquid pressure drop to the all-gas pres- 


sure drop is the square of the two-phase flow modulus, ca or 


7 
a 
XY _ AF 
= ABP; ei 
Figure 23 shows the relationship between Y and the correlation param- 
eter @ for each of the various flow regimes, as determined experiment- 
ally by Lockhart and Martinelli. If @, is taken from the curve the 


two-phase pressure drop is 


eo 
Nees = Bb, AR. (45) 
Similar ly, if , is taken from the curve the two-phase pressure drop is 


2 
APre = b, AF, (46) 


It should be noted that Lockhart and Martinelli state that 
more data are needed to establish the validity of their correlation at 
very high and very low values of the flow modulus x (nearly all- 


liquid and nearly all-gas flow, respectively). 
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2. Chenoweth-Martin Two-Phase Flow Correlation 

Chenoweth and Martin developed an "improved correlation for 
two-phase pressure drop in horizontal pipes’ which is especially valid 
for pressures up to 100 psia. The correlations can be used for any 
two-phase mixture as long as the flow is turbulent and it predicts 
single-phase values when the flow is all liquid or all gas. The two 
imput parameters for the use of the correlation are the liquid volume 
fraction (LVF) and the ratio of a fictitious all-gas pressure drop to 
a fictitious all-liquid pressure drop. 

The liquid volume fraction is calculated from the flow rates 


and densities of the two-phases using the equation 


; : Me 4 
LVF = Volume flow rate of liquid _ /°.. ‘ 


Total volume flow rate ef, + 6/0. 1+ tale 


The ratio of the fictitious pressure drops is determined by computing 


the value of Pr De"/ pc yp * 


where 


Dx Pee + Lk 


(48) 


and 
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The subscripts L* and G* denote the fictitious all-liquid and all-gas 
states, respectively. The K's are loss coefficients for valves and 
fittings. The fictitious friction factor is determined in the normal 
manner using an artificial Reynolds number based on the total flow rate 
of liquid and gas and the physical properties of the phase in question. 
is determined using the Reynolds number 

Re,» = ed (50) 
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The fictitious all-liquid pressure drop is then computed as 


For example, fi 


2 
> it 
Ars = phx 5 z (51) 


where V, .. is a fictitious velocity determined from the equation 
m1 + me 
PLA (52) 


Figure 24 can now be used to predict the two-phase pressure drop, AP-p. 


x = 


i 


In contrast to the Lockhart-Martinelli correlation, the 
Chenoweth-Martin correlation covers only all-turbulent two-phase 
mixtures but ties into single-phase values at either end. It is 


also valid over a wider range of presSures. 


B. EQUIPMENT DESIGN OBJECTIVES AND DESCRIPTION 

A test rig was developed for the experimental analysis of two-phase 
flows, in particular the mixing process which immediately follows 
injection of water into the air. Initially, fully established two- 
phase flows will be investigated by directly measuring the wall fric- 
tion force and the overall pressure drop. Subsequent investigation 
will concentrate on determination of the pressure drop due to the mixing 
process in a developing flow. The direct measurement of the wall fric- 
tion force will allow separation of the pressure drop due to mixing from 
the overall pressure drop. General views of the apparatus are shown 
in Figures 25 and 26. A flow diagram for the initial tests is shown 
in Figure 27. 

Air enters the system through a 2-inch pipe from a plenum fed by a 
reciprocating compressor capable of supplying air at pressures up to 
150 psig and flow rates up to /00 scfm. Water is supplied to the sys- 
tem through a 3-inch pipe from a centrifugal pump capable of supplying 
up to 500 gpm at 80 psig. The mass flow rates of the air and water are 
accurately determined using thin-plate orifice meters. The air and 
water flows are united by joining the two pipes into a 3-inch pipe. 

The two-phase mixture flows together for about 28 feet becoming fully 
established, then flows through a "floating'' test section 10 feet long 
in which the wall shear force is measured. The downstream end of the 
test section exhausts to the atmosphere where the water is collected 
and returned to the water pump. 

1. Flow Orifice Details 

The flow measurement orifices shown in Fig. 28 were designed 
according to standard ASME specifications [Ref. 10] . Figures 29A 
and 29B give the dimensions of the orifices and the placement of the 
pressure taps with respect to the pipe. Each orifice plate was 
machined from type 304 stainless steel plate for which a coefficient 


of expansion curve is available in Ref. 10. Each of the flanges con- 
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taining an orifice was machined with a circular groove centered with 
respect to the pipe axis to insure accurate alignment of the orifice 
with the pipe axis. 

2. witial Mixwme of Air and Water 

The mixing of the air and water is accomplished by merely 
joining the two pipes at a 45-degree angle (see Fig. 30). The mixing 
process occurs solely due to the turbulent nature of the flow; if the 
liquid volume fraction is sufficiently small a mist-flow will result. 
It is the mist-flow regime (liquid volume fractions below about 15 
percent) which is directly applicable to the water-augmented turbofan 
engine described in Section I. However, the apparatus is designed to 
handle as wide a range of water-to-air mass flow ratios as permitted 
by the available water pump and air compressor capacities. 

3. Test Section Details 

The test section consists of a 10-foot long, 3.106-inch steel 
pipe which is supported two feet from each end. Each support has four, 
one-eighth inch diameter stainless steel rods, one vertically above, one 
vertically below, and one horizontally on each side as shown in Fig. 31. 
The rods are mounted on bearings which allow the test section to "float" 
Veneitudinally. 

The wall shear force measurement is made through two cantilever 
flexures mounted horizontally, one on each side of the test section 
(see Fig. 32). Each flexure has four strain gages mounted to detect 
changes in bending moment caused by the force exer”~ed by the test sec- 
tion. Figure 33 shows the wiring diagram for each flexure. 

The labyrinth seal shown in Fig. 34 was designed to provide a 
pressure seal at the inlet to the test section which does not restrict 
longitudinal motion. The dimensional details of the labyrinth are 
given in Fig. 35. Experimental determination of the leakage rate is 
discussed in Section C below. 

The pressure at the inlet to the test section is measured to 
allow the determination of the pressure drop through the test section. 
Additionally the test section inlet pressure is required in order to 
calculate the force on the labyrinth face at the test section inlet. 
This pressure force must be subtracted from the total force indicated 


on the flexures to obtain the true wall shear force. 
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C. CALIBRATION PROCEDURE 

It is necessary to determine the amount of air leakage which occurs 
at the labyrinth seal over the range of expected internal to external 
pressure ratios. This was to have been accomplished by sealing the 
outlet of the test section, varying the internal pressure and measuring 
the mass flow rate of air leaving through the labyrinth seal. However, 
excessive leakage in the air supply line prevented the accurate meas- 
urement of the leakage rate through the labyrinth seal and insufficient 
time remained to complete the necessary repairs. 

The flexures on which the strain gages are mounted require calibra- 
tion before and after each series of test runs. The calibration proced- 
ure for the flexures involves placing a known force on the test section 
using a system of weights and a pulley as shown in Fig. 36. A typical 
calibration curve of force versus strain gage output is presented in 
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D. DATA COLLECTION AND REDUCTION PROCEDURE 

The important variables for the operation of the two-phase flow 
test rig are the mass flow rates of air and water, the force exerted 
on the test section due to wall friction and labyrinth pressure effects, 
and the pressure loss in the test section. To determine the mass flow 
rates it is necessary to record the static pressure upstream of the air 
orifice, the pressure drop across both orifices for both the flange 
taps and the D and half-D taps, respectively, the temperature of each 
fluid and the atmospheric temperature at each manometer. 

The calculation of the two-phase pressure drop predictions of both 
Lockhart and Martinelli and Chenoweth and Martin require that the inlet 
Static pressure in the test section and the temperature of the air-water 
mixture at the outlet of the test section be recorded. Atmospheric 
pressure must be recorded for use in mass flow rate calculations. The 
temperature is recorded at the outlet of the test section because the 
temperature probe will disturb the flow if it is installed at the inlet 
to the test section. It should be mentioned at this point that the 
temperature measured in a two-phase mixture can be subject to discussion 
because it is not known what effect the water droplets impinging on the 


probe have on the indicated temperature. The accuracy of the pressure 
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drop predictions depends heavily upon an accurate temperature of 
the flow. 

A computer program converts the raw experimental data into mass 
flow rates and two-phase pressure drop predictions. The input and 
output variables used for the computer program are defined in Appendix 
III and the computer program itself along with a sample output is pre- 
sented in Appendix IV. 

The pressure drop prediction of Chenoweth and Martin is not con- 
pletely carried out by the computer program because of the difficulty 
in obtaining an analytical expression for the curves of the correlation 
in Fig. 24. The program does supply the variables necessary to enter 
Fig. 24. The value of AP-p /AP,x thus found is multiplied by the 
fictitious all-liquid pressure drop, AP x , which is also computed 


in the program, to obtain the predicted two-phase pressure drop. 


E. CURRENT STATUS AND SUGGESTIONS FOR FURTHER WORK 
The work reported herein concerning the development of the two- 
phase flow test rig suffered several setbacks due to difficulty in 
obtaining parts. Because of these delays the installation of the 
water return system has not been completed. 
Several suggestions for further work have become apparent during 
the course of the development of the two-phase flow test rig. 
1. Labyrinth Air leakage Determination 
It will be necessary to determine the leakage rate of air 
through the labyrinth seal prior to taking experimental data. It 
should be noted that the possibility exists that the present metering 
system may not be adequate for measuring minimum flow rates of both 
air and water. 
2. Labyrinth Water Leakage Determination 
During operation with two-phase or all-liquid flow some water 
will probably leak out the labyrinth seal and provision must be made 
to measure this leakage. 
3. Computerization of Friction Factor Curve 
The friction factors used in calculating pressure drops in the 
computer program solutions for predicted two-phase pressure drop were 
assumed to be constant values. This is true only in the region of 


fully turbulent flow. The value of this friction factor used for high 
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Reynolds numbers should first be checked by simple pressure drop tests. 
The computer program must then be improved by the computerization of 
the friction factor curve so that the laminar flow and transition re- 
gions will be included. 
4. Computerization of Chenoweth-Martin Correlation Curve 
It would be helpful if an analytical expression could be found 
for the Chenoweth and Martin correlation curves of Fig. 24. 
+S. Hvaluationjer Predicted) Pressure soup 
The knowledge gained from the continuation of the work reported 
herein should be used to evaluate the pressure drop correlations of 
Lockhart and Martinelli and Chenoweth and Martin. Also, a prediction 
for the two-phase wall shear has very recently been proposed by Muench 
and Ford [Ree. 11] and should be evaluated. The analysis of the two- 
phase flow calculations in the water-augmented turbofan should be re- 
examined in the light of experimental results. 
6. Pressure Drop due to Two-Phase Mixing 
The present test rig is set up to obtain fully developed two- 
phase flow in the test section. However if the water is injected im- 
mediately upstream from the test section it should be possible to 
separate the pressure drop into that caused by the mixing of the two 
phases and that caused by wall friction. This knowledge would be use- 


ful in further analysis of injection systems. 
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FIGURE 2 .THRUST RATIO VERSUS FAN PRESSURE RATIO 
FOR VARIOUS WATER INJECTION VELOCITIES 
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FIGURE 3.THRUST RATIO VERSUS FAN PRESSURE RATIO 
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FIGURE 4 . THRUST RATIO VERSUS BYPASS 
RATIO FOR VARIOUS INJECTION VELOCITIES 
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FIGURE 5 . THRUST RATIO VERSUS BYPASS RATIO 
FOR VARIOUS WATER-TO-GAS GENERATOR AIR RATIOS 
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FIGURE 6 . THRUST RATIO VERSUS WATER-TO-GAS 
GENERATOR AIR RATIO FOR VARIOUS INJECTION VELOCITIES 
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FIGURE 7.EFFICIENCY RATIO VERSUS WATER-TO-GAS 
GENERATOR AIR RATIO FOR VARIOUS INJECTION VELOCITIES 
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RATIO FOR VARIOUS WATER INJECTION VELOCITIES 
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FIGURE |O .EFFIENCY RATIO VERSUS FAN DUCT VELOCITY 
RATIO FOR VARIOUS WATER INJECTION VELOCITIES 
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FIGURE I! .THRUST RATIO VERSUS FAN DUCT TEMPERATURE 
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FIGURE 13.THRUST RATIO VERSUS FANNOZZLE VELOCITY 
RATIO FOR VARIOUS WATER INJECTION VELOCITIES 
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FIGURE 21 EFFICIENCY RATIO VERSUS FAN MIXING DUCT 


PRESSURE RATIO FOR VARIOUS WATER INJECTION 
ji VELOCITES 


61 


o 
O 


vr 
O 


ee ee 8 he ee Dh eh 


eal 


O 


ps ne 
& 30 
a 
5 — 
ee ee 
20 —o S/S 
: re-ltine 
uJ 
: ; = = 
¥ a] 


BF 
3 
i 


O 


O9C 1.00 Ltd 


Fan Mixing Duct Pressure Ratio — Pag / Pa7 


Oo 
3 
fo 


> 
O 
oe | 


FIGURE 22.EFFCIENCY RATIO VERSUS FAN MIXING DUCT 
PRESSURE RATIO FOR VARIOUS WATER-TO-GAS GENERATOR 
AIR RATIOS 7 


62 





@ ‘19J9WOUOY UOYWO/GI0D MO]4 8SOYd — OML 


63 


| Two-Phase Flow Modulus , xX 


FIGURE 23.TWO PHASE FLOW CORRELATION OF LOCKHART—MARTINELLI (REF. 5] 
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FIGURE 31 


TEST SECTION SUPPORTS 


71 





ta 


FIGURE 32 


FLEXURES 








Fine Balance 


Coarse Balance 


Exitation (-) 


A. Switching And Balancing Unit 


B. Daytronic Strain Gage Digital Indicator 
Mode! 700 


C. Compression Strain Gages 
T. Tension Strain Gages 


Strain Gages: Budd 12042 Foil Gages 
Bonding Agent. Eastman 910 


Waterproofing: Dow Corning RTV 


FIGURE 33 
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FLEXURE CALIBRATION SET-UP 
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APPENDIX I 


TURBOFAN PROGRAM NOMENCLATURE 


Input Parameters 
Card le 


Gas generator diffuser efficiency 
Compressor efficiency 

Burner efficiency 

Turbine efficiency 

Gas generator nozzle efficiency 
Fan duct diffuser efficiency 

Fan efficiency 

Fan nozzle efficiency 

Water injection pump efficiency 


Ambient density (slug/£5~ 


Ambient pressure (lb/ft ) 
Ambient temperature (degrees Rankine) 
Heating value of fuel (BIU/1bm) 

Reference air temperature (492 deg. R) 


) 


Air injection mach number (at fan duct entrance) 
Sea water temperature (deg. R) 

Burner total pressure ratio 

Fan duct total pressure ratio (dry fan) 


Number of vessel! speeds 

Number of compressor pressure ratios 
Number of fan pressure ratios 

Number of turbine inlet temperatures 
Number of bypass ratios 

Number of water-to-gas generator air ratios 
Number of mixing duct pressire ratios 
Number of mixing duct velocity ratios 
Number of mixing duct temperature racios 
Number of fan nozzle velocify ratios 
Number of fan nozzle: temperature ratios 
Number of water injection velocities 


Cool temperature specific heat ratio 
Hot temperature specific heat ratio 
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Gard 6: 
VOCL) 


Gard /: 
PTR21(J) 
Gand, 8; 
PTR76(M) 
Gard? 9 ; 
TI3(K) 
Gard 10; 
BR(L) 
Gard 11; 
WGR(N) 
Card Iz; 
PSReae 
Card 13: 
VR8WACIK) 
Card 14; 
TSR8WACIL) 
Canad 15; 
VROWACIM) 
Card 16: 
TSROWACIN) 
Gard 17; 
VPW(NN) 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


of 


OF 


of 


of 


of 


OF 


OF 


Of 


fee 


OF 


of 


of 


vessel speed (knots) 


compressor pressure ratio 


fan pressure ratio 


turbine inlet temperature (deg. R) 


bypass ratio 


water-to-gas generator air ratio 


mixing duct static pressure ratio 


mixing duct velocity ratio 


mixing duct temperature ratio 


fan nozzle velocity ratio 


fan nozzle temperature ratio 


water injection velocity (knots) 


OUTPUT PARAMETERS NOT PREVIOUSLY DEFINED 


PWORK 
DHT43 
STHRUS 
ER 
ETAPRO 
ERATIO 


Water injection pump work (ft-lb/slug) 
Turbine work (ft-lb/slug) 


Specific thrust (1b thrust/1bm/sec) 


Ratio of water-augmented-to-dry specific thrust 


Propulsive efficiency 


Ratio of water-augmented-to-dry propulsive efficiency 


io 


SFC 
VA/ 
TS8 
VA8 
og 
VAY 


Specific fuel consumption (1b fuel/1lb thrust-hr) 
Air velocity at fan duct entrance (ft/sec) 

Air static temperature at mixing duct exit (deg. R) 
Air velocity at mixing duct exit (ft/sec) 

Air static temperature at fan nozzle exit (deg. R) 


Air velocity at fan nozzle exit (ft/sec) 
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APPENDIX IIT 


DATA REDUCTICN PROGRAM NOMENCLATURE 


INPUT PARAMETERS 


Ma Number of interpolation pairs for subroutine SPLIN1 
R_ _ (I) Reynolds number values 
K_= (I) Discharge coefficient values corresponding to R__ (1) 
pubscrip@s Heme oa. kh toe kK Ct). 
FA Flange pressure taps, air 
DA D and % D pressure taps, air 
FW Pillange pressure taps, water 
DW D and % D pressure taps, water 
Example: RDA(I) Reynolds number for D and 4% D taps for air flow 
hal Number of interpolation pairs for subroutine SPLIN1 
Xo 3GE) Lockhart-Martinelli two-phase flow modulus values, % 
F___-_ ~ (I) Lockhart-Martinelli correlation parameter values, Dd 
Sugeeripts for M , X (i,m £1): 
G Lockhart-Martinelli gas correlation 
L Lockhart-Martinelli liquid correlation 
Te Turbulent-turbulent flow 
TV Turbulent-viscous flow 
VT Viscous-turbulent flow 
VV Viscous-viscous flow 
Example: FGIV(I) Correlation parameter, d , for gas-phase 
turbulent-viscous flow 
DIA Inside diameter of air pipe at. orifice (in) 
D2A Diameter of airvonlmieomii) 
D1W Inside diameter of water pipe (in) 
D2W Diameter Of Wale. Orliree cin) 
HMFA Manometer fluid differential, flange taps, air (in H,0) 
HMDA Manometer fluid differential, D and 4 D taps, air (in H,9) 
HMF W Manometer fluid differential, flange taps, water (in Hg) 
HMDW Manometer fluid differential, D and 4 D taps, water (in Hg) 
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TLFA 
TLFW 
TAFA 
TAF W 
El 
Pa 
TIFA 
TIF W 
DT 
LT 


Air flow temperature at orifice (deg. F) 

Water flow temperature at orifice (deg. F) 

Air temperature at air orifice manometer (deg. F) 
Air temperature at water orifice manometer (deg. F) 
Air static pressure upstream of orifice ae 
Static pressure at test section inlet one 

Air temperature at tesx section outlet (deg. F) 
Water temperature at test section outlet (deg. F) 
Inside dieeter of test section (in) 


Length of test section (ft) 


OUTPUT PARAMETERS NOT PREVIOUSLY DEFINED 


BETAA 
BETAW 


RFAL 
RDAIL 
RFWIL 
RDWI 
KFAIL 
KDAT 
KF WI 
KDW1L 
FRFA 
FRDA 
FREW 
FRDW 


RHOWT 
RHOAT 


RLM 
RGM 
VL 
VG 


PDROPL 
PDROPG 
eee Di 


Ratio of orifice-to-pipe inside diameter, air pipe 
Ratio of orifice-to-pipe inside diameter, water pipe 
Air Reynolds number based on flange taps 

Air Reynolds number based on D and 4 D taps 

Water Reynolds number based on flange taps 

Water Reynolds number based on D and ¥ D taps 

Air discharge coefficient based on flange taps 

Air discharge coefficient based on D and % D taps 
Water discharge coefficient based on flange taps 
Water discharge coefficient based on D and % D taps 
Air mass flow rate based on flange taps (l1lbm/sec) 
Air mass flow rate based on D and 4 D taps (l1bm/sec) 
Water mass flow rate based on flange taps (lbm/sec) 
Water mass flow rate based on D BA 45 D taps (lbm/sec) 
Inside diameter of test section (ft) 

Water density in test section @ienfee) 

Air density in test section (lpmye te) 

Liquid Reynolds number (Lockhart-Martinelli) 

Gas Reynolds number (Lockhart-Martinelli) 

Liquid velocity (ft/sec) 

Gas velocity (ft/sec) 

Fictitious liquid pressure drop (1b/£t"/10£t) 
Fictitious gas pressure drop (lice /106t ) 


Lockhgrt-Martinelli two-phase pressure drop prediction 


Cloyne Ort) 


23 


RLC Liquid Reynolds number (Chenoweth-Martin) 


RGC Gas Reynolds number (Chenoweth-Martin) 

PDLC Dimensionless group, p, (all-liquid flow) 

PDGC Dimensionless group, We (all-gas flow) 

CORPAR Chenoweth-Martin correlation parameter 

VFR Volumetric flow ratio (liquid-to-gas) 

PDRLC Fictitious all-liquid pressure drop (1b/£t7/10£t) 


LVOLFR Liquid volume fraction 


DATA INPUT CARD ORDER 


Card 1: MFA, MDA, MFW, MDW 
Card 2: RFA(I), I = 1, MPA 
Card 3:  KFA(I), = 1, MFA 

Card 4: RDA(I), I = 1, MDA 
Card 5: akKDACL), a=" 
Card 6: RFW(I), I = 1, MFW 
Card 7: KFW(I), I = 1, MEW 
Card 8: RDW(I), I = 1, MDW 
Card 9: KDW(IT), 1 = 1, MDW 
Card 10: DIA, D2A, D1W, D2W 
Card 11: HMPA HMDA, TIPFA, TAPA, Pl 
Card 12: HMFW, HMDW, TIFW, TAFW 
Card 13: PT, TTFW, TIFA 

Card 14: DT, LT 

Card 15: MG, ML 

Card 16: XG(I), I = 1, MG 

Card 17: XL(I), I= 1, ML 

Card 18: FGTT(I), I = 1, MG 
Card 19: FGTV(I), I = 1, MG 
Card 20: FGVT(I), I = 1, MG 
Card 21: FGVV(I), I = 1, MG 
Card 22: FLTT(1), I = 1, ML 
Card 23: FLTV(I), I = 1, ML 
Card. 24: FLVT(I), I = 1, ML 
Gard 25: “PEVV( 1) esi 
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WATER AND AIR MASS FLOW RATE CETERMINATION 


SQUARE ECGE CONCENTRIC ORIFICE WITh FLANGE TAPS AND DEHALF-D TAPS 
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STANDARD ASME METHCOS USED 


ORLEICE MATERIAL 


TYPE 304 STAINLESS STEEL 
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CALCULATION OF LOCKHART-MARTINELLI AND CHENOWETH-MARTIN 


TWO-PHASE PRESSURE CROP PREDICTION 
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